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SYNCHRONOUS  DEDICATED  PROCEDURES 


Activation 

1.  Power  on  equipment. 

2.  When  connected,  identify  your  terminal  to  the  network  by  entering 

0$$S0N  XXXXXX  (enter  your  slgn-on  identification  code  for  XXXXXX) 

(Examp:  e:  t>$$SON  AA1234)  1 1 

3.  The  system  will  respond  with: 

UNIVAC's  Tel con  Network-Level  XXX  DCPID-XXX 
SESSION  PATH  OPEN 

4.  Identify  your  site  to  the  EROLS  system  by  transmitting  the  following: 
C>SGNONS/TERMINAL  ID 

(Example:  DSGNONS/ABCDE)  £>  5£rA/0  KSjflTlS  (j>£> 

5.  Thf  system  will  respond  with: 

*MSG  ONI  SIGN-ON  ACCEPTED 

6.  Enter  your  terminal  ID  as  in  the  past.  Entering  of  the  terminal  ID  at  this 
point  will  be  eliminated  in  the  near  future. 

7.  DROLS  commands  remain  the  same. 

Termination 

1.  Enter  @TERM@  and  transmit. 

2.  The  system  will  respond  with: 

THIS  TERMINAL  HAS  BEEN  TERMINATED 
*MSG  D07  PLEASE  SIGN  OFF  TERMINAL 

3.  Disconnect  (sign  off)  from  the  network  by  transmitting: 

>$$SOFF 

4.  The  system  will  respond  with: 

SESSION  PATH  CLOSED 

5.  Power  off  equipment. 


L.' 


P 


s 


\  ■ 

i 

[  ' 

n 


r 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wbon  Data  Bntorod) 


REPORT  DOCUMENT ATIOM  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  report  number 

TECHNICAL  REPORT  NO.  84 

2,  GOVT  ACCESSION  NO. 

J.  RECIPIENT’S  CATALOG  NUMBER 

M 

4.  TITLE  (and  Subtitle) 

DIOPHANTINE  EQUATIONS  AND  THE  CUBIC  AMBIGUITY 
RESOLUTION  ALGORITHM 

S.  TYPE  OF  REPORT  *  PERIOD  COVERED 

8.  PERFORMING  GRG.  REPOilT  NUMBER 

7 

7.  author/*; 

JAMES  T.  CATER,  JR. 

8.  CONTRACT  OR  GRANT  NUMBER/*; 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Data  Sciences  Division 

National  Range  Operations  Directorate 

White  Sands  Missile  Range,  New  Mexico  88002 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  8  WORK  UNIT  NUMBERS 

1 1,  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Mathematical  Services  Branch  (STEWS-NR-AM) 

12.  REPORT  DATE 

OCTOBER  1982 

Data  Sciences  Division,  National  Range  Operations 
White  Sands  Missile  Range,  NM  88002 

13.  NUMBER  OF  PAGES 

25 

i 

U.  MONITORING  AGENCY  NAME  *  ADDRESS/!/  <11  Heron  t  from  Controlling  OlHeo) 

15.  SECURITY  CLASS,  (ot  Mo  roport) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/ DOWNGRADING 

N/A 

1*.  DISTRIBUTION  STATEMENT  (ot  Mo  Roport) 

l 

■ 

DISTRIBUTION  OF  THIS  DOCUMENT  IS  UNLIMITED. 

5 

'7.  DISTRIBUTION  STATEMENT  (of  thm  obotrect  ontotod  In  Sloe*  20,  II  dllloront  Item  Roport) 

U.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  tmrerao  aide  it  neceeeery  mtd  identity  by  block  number) 

Resolution 

Integer 

Diophantine 

20.  ABSTRACT  (Contlnuo  on  rotroroo  mi  do  II  noc  oooory  and  Idontlty  by  block  numb#r) 

i 

' )  The  ambiguity  resolution  problem  for  the  TOAME  four  phase  interferometer 
system  is  presented  in  terms  of  the  integer  solutions  of  systems  of  linear 
equations 

DD  ,  jSnM73  1473  EDITION  OF  I  NOV  «5  IS  OBSOLETE  UNCLASSIFIED 


t 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Hhit  Doto  Eutorod) 


The  ambiguity  resolution  problem  for  the  TOAME  four  phase  interferometer 
system  is  presented  in  terms  of  the  integer  solutions  of  systems  of  linear 
equations . 


TABLE  OF  CONTENTS 


INTRODUCTION  . 

A  BRIEF  REVIEW  . 

THE  OIOPHANTINE  APPROACH 
THE  CUBIC  ALGORITHM.  .  . 


INTRODUCTION 


This  report  deals  with  data  from  the  TOAME/PDM  systems  at  11-si  te  and  A-site, 
A  Oiophantine  equation  approach  to  ambiguity  resolution  for  four  phases  is  pre¬ 
sented.  Also  the  currently  used  Cubic  algorithm  is  shown  to  be  a  variant  of 
this  approach.  Both  algorithms  are  exercised  against  difficult  data  recorded 
at  the  sites  from  a  high  performance  vehicle.  Plots  of  cosine  versus  time  are 
included.  They  are  self-explanatory.  The  algorithm  presented  here  performs 
better  than  the  Cubic  algorithm  on  this  data. 

A  BRIEF  REVIEW 

Some  familiarity  with  the  problem  and  the  nature  of  the  data  is  assured.  At 
each  site  there  are  two  orthogonal  base  lines.  Each  base  line  contains  four 
pairs  of  antennas  with  a  phase  meter  for  each  pair.  The  distances  spanned  by 
each  pair  are  ideally  given  by  the  vector  bT  *  (43,  42,  36,  280)  where  1  <  A/2 
(slightly  less  than  one-half  tne  wave  length  of  the  received  signal).  Unfor¬ 
tunately  the  actual  spacings  of  the  antennas  do  not  conform  too  well  to  these 
proportions  thus  introducing  modeling  errors.  The  quantities  to  be  ultimately 
found  are  the  differences  in  range  of  the  target  between  the  two  antennas  or 
each  pair.  These  range  differences  are  measured  in  cycles  of  the  operating 
frequency.  Designate  these  range  differences  by  the  vector  dT  =  ( d x ,  d2,  d3, 
d4).  The  scaled  output  of  the  phase  meters  p.,  j  =  1,  .  .  ,  ,  4  is  related  to 
the  components  of  dT  by  the  equations  d.  =  n.  +  p.,  j  *  1 . 4  where  the 

J  J  J 

n.  are  positive  or  negative  integers  and  0  <  p.  <  1.  The  various  algorithms 

J  >' 

attempt  to  recover  tne  n.  using  the  p.  as  inputs. 

vJ'  J 
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THE  DIOPHANTINE  APPROACH 

Again  let  b^  *  (43,  42,  36,  285)  and  d  *  n  +  p  where  n  and  p  are  vectors  con¬ 
sisting  of  the  components  n.  and  p.  respectively.  The  planar  wave  front 

J  0 

assumption  implies  that 
d  »  yb 

or  equivalently 

n  +  p  *  yb,  where  y  is  a  scalar. 

Now  there  is  an  abundance  of  3  by  4  integer  matrices  A  such  that 
Ab  =»  0 

Hence 

A  (n  +  p)  3  yAb  3  0 

or 

An  =  -Ap 

The  left  side  of  the  last  equation,  An,  is  clearly  an  integer  vector  so  it  is 

expected  that  the  right  side,  -Ap,  is  also  an  integer  vector.  If  this  were  true 

then  the  system  of  Diophantine  equations 

An  =  -Ap  =  m 

could  by  solved  for  the  vector  n.  However  in  reality  -Ap,  when  computed,  is 

never  an  integer  vector  regardless  of  what  suitable  matrix  A  is  used.  The  rea¬ 

sons  for  this  include  the  facts  that  the  antennas  are  not  quite  correctly 


positioned,  che  components  of  p  are  stored  in  finite  data  registers,  and  the 
planar  wavefront  assumption  is  not  exact.  Lump  any  errors  or  discrepancies 
in  the  vectors  b  and  P  into  the  error  vector  e.  Mow  rewrite  the  equation  as 
An  a  -A(p+e)  -  m+Ae 

where  the  error  vector  e  is  unknown.  About  the  only  thing  is  to  attempt  to 
obtain  the  integer  vector  m  by  means  of  the  nearest  integer  operator.  Denote 
the  nearest  integer  of  x  by  ANINT(x)  in  accord  with  FORTRAN.  When  the  argument 
is  a  vector  the  intent  is  a  component-wise  operation.  This  leads  to  the 
Oiophantine  system 
An  3  -ANINT(A  (p+e))  =  ANINT  (m+Ae) 

Even  though  ANINT  is  a  nonlinear  operator,  it  is  true  that 
ANINT  (integer  +  real  number)  3  integer  +  ANINT  (real  number) 

Thus 

An  3  ~ANINT(A  (p+e))  *  ANINT  (  Ae)  +  m 

If  any  component  of  ANINT  (Ae)  is  nonzero  then  an  error  is  present.  Worse 
yet  this  error  is  amplified  in  the  solution  of  the  system  of  equations.  These 
errors  cause  spikes  and  step  discontinuities  in  the  cosine  plots.  Some  of 
the  enclosed  cosine  plots  illustrate  this  phenomenon. 

There  probably  is  no  certain  way  of  squashing  ANINT(Ae)for  all  conditions,  but 
we  can  try  to  minimize  the  magnitudes  of  the  components  of  Ae.  One  way  to  get 
a  handle  on  this  problem  and  also  state  it  more  elegantly  is  to  say  we  want  to 
minimize  the  max-norm  (alias  °°-norm)  of  Ae.  The  max-norm  of  a  vector  is  the 
absolute  value  of  the  component  having  the  largest  magnitude.  So  it  becomes 
apparent  that  we  want  a  matrix  A  whose  components  have  small  magnitudes  or  we 
would  like  to  minimize  the  induced  norm  of  A.  We  may  not  actually  minimize  this 
induced  norm  (max  l  |a.J)  but  we  try  to  make  it  small  and  hope  for  the  best. 

1  j  J 


Start  with  trial  matrix  AQ  such  that  AQb 
42  -43  0  0 

Aq  =  0  6  -7  0 

0  0  8  -1 


For  instance  let 


A  direct  attack  cn  the  system  A Qn  =  -ANINT(Agp)  encounters  awkward  difficulties, 
It  is  prudent  to  factor  Ag  in  the  following  way: 
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Without  getting 

involved 

in 

details , 

the  equation 

V  = 

-A  p  can 
o  v 

be  lef 

multiplied  by  appropriate  matrices  resulting  in  A^n  «  -A-j  p  where 
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Notice  that  still  A^b  =  0.  Furthermore  the  equation  A-|  n  =  -AN  I  NT  ( p) 
does  have  an  integer  solution.  However  the  top  row  of  A^  contains  components 
of  rather  large  magnitudes.  In  order  to  solve  this  Diophantine  system  we  will 
disguise  somewhat  the  determination  of  the  extra  parameter  caused  by  having 
3  equations  4n  4  unknowns.  This  is  done  by  an  indirect  use  of  the  fictit-'ons 
"very  coarse"  measurement.  Let  B.|  designate  the  righthand  4  by  4  matrix. 


» 


That  is  to  say 


-1  0 


8  -1 


Now  let  m  =  -ANINT(B^p),  Next  solve  the  system  B^n  =  m.  The  solution  is 

"  i  _  “i 

given  by  n  =  3^  m.  Since  det(B^)  =  1  it  is  guaranteed  that  5^  is  an 
integer  matrix  and  the  solution  vector  n  has  integer  components.  It  is  only 
the  last  component  of  n  that  is  of  interest.  This  is  given  by 
n^  =*  -8  (m.|  +  m^)  -  m3  +  288  m^.  The  graph  entitled  'TEST'  shows  the  plot 
of  (n^  +  p^)/(base  in  feet)  or  cos  against  time.  This  is  a  real  disaster  and 
we  suspect  the  large  elements  in  the  first  row  of  or  equivalently  B^ . 

Every  detail  is  not  being  explained  here  but  we  are  going  to  use  elementary 
row  operations  to  transform  the  first  three  rows  of  B-j .  We  are  interested  in 
only  those  row  operations  which  do  not  change  the  magnitude  of  the  determinant. 
This  means  we  can  replace  the  i—*1  row  by  minus  the  i—  row  and  we  can  replace 
the  i™  row  by  i—  row  =  i—  row  +  integer  *  j—  row.  We  apply  these  operations 
with  the  objective  of  eliminating  large  magnitudes  among  the  elements  of  B^ . 
Specifically  make  the  following  simultaneous  replacements:  1  —  row  = 
l5^  row  +  7  *  2—  row  +  3—  row  and  2—  row  =  -1—  row  -6*2—  row  -3—  row. 

This  is  tantamount  to  left  multiplying  B-^  by  ^  ^  ^ 

-1  -6  -1  0 

0  0  10 

0  0  0  1 
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Call  the  resulting  matrix  B2  where 


0  6  1-1 

-6  1-2  1 

0  0  8  -1 

1-1  0  0 


Since  det(B2)  =1  it  is  assured  that  B2  is  also  an  integer  matrix.  Now  let 
m  =  -ANINT(B2P).  Solve  B2n  =  m  for  the  integer  vector  n  =*  B2  in.  The  last 
component  of  n  is  now  given  by 

n4  =  40m-j  +  48m2  +  7m3  +  288m4 


The  graphs  of  the  cosine  using  this  algorithm  are  entitled  1 DIOPH ' .  It 
appears  to  work  quite  well  on  the  relatively  difficult  data  used  here.  Each 
graph  of  ' DIOPH '  is  followed  by  the  corresponding  graph  of  the  Cubic 
algorithm.  One  notable  fact  is  that  where  the  Cubic  plot  is  smooth  the 
quantity  n4  +  p4  is  bit  by  bit  the  same  in  both  algorithms. 

The  Cubic  Algorithm 

The  Cubic  algorithm  is  currently  in  use.  In  order  to  describe  it  first  define 
the  following  two  composite  inputs. 

Pq  =  Pi  -  P2  which  corresponds  to  a  fictitious  base  lire  of  length  unity.  This  is 
called  the  very  coarse  input. 

Pc  =  P2  -  P3  which  corresponds  to  a  fictitious  base  line  of  length  6.  This  is 
called  the  coarse  input. 

Now  the  algorithm  can  be  described  as  a  three  stage  procedure. 
nc  =  ANINT  [6Po  -  Pc3 
n3  =  ANINT  [6(nc  +  Pc)  -  ?3] 
n4  ■  ANINT  [8 ( n 3  +  P3)  -  P4] 
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The  quantities  n3  and  n4  have  exactly  the  same  meanings  as  in  the  Oiophantine 
solution  and  n4  is  the  quantity  of  interest.  At  this  point  recall  that 
ANINT  (integer  +  real  number)  =  integer  +  ANINT  (real  number) 

Apply  this  to  the  expression  for  n^. 
n4  =  8n3  +  ANINT(8p3  -  p4) 

Substitute  for  n3. 

n4  =  48nc  +  8  *  ANINT(6pc  -  p3)  +  ANINT(8p3  -  p4) 

Continuing. 

n4  =  48  *  ANINT(6n0  -  pc)  +  8  *  ANINT(6pc  -  p3)  +  ANINT(8p3  -  P4) 

Substitute  for  pc< 

n4  =  48  *  ANINT(6dq  -  P2  +  P3)  +  8  *  ANINT(5p2  -  7P3)  +  ANINT(8p3  -  P4) 

For  future  convenience  rearrange  signs. 

n4  =  -48  *  ANINT(-6p0  +  p2  -  p3)  -  8  *  ANINT(-6p2  +  7p3)  -  ANINT(-8P3  +  P4) 


In  the  description  of  a  Diophantine  eq<  ^tion  approach  we  started  with 


T 

n 


A0n  =  -AqP 


( ^i  >  n2 » ,  n4 )  >  p  =  ( p^  >  p2 « p3 p  4 ) »  Aq  - 


42  -43 

0  6 

0  0 


0  0 

-7  0 

8  -1 


Augment  this  system  with  the  very  coarse  measurement  pg  =  p-j  -  p2.*  Now  we  have 
nT  =  (nQ,n1 ,n2,n3,n4)  and  pT  =  (pg,p1 ,p2,p3>p4) 

The  measurement  Pg  corresponds  to  a  fictitious  base  line  bg  where  bQ  «  1 . 

Since  1 <A/2  it  follows  that  always  nQ  =  0.  Furthermore  make  the  redefinitions 


b  =  (bg,b1  ,b2,b3,b4)  and  Arj  = 


43 

0 

0 

0 


-1 

42 

0 

0 


0 

-43 

6 

0 


0 

0 

-7 
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0 

0 

0 

-1 
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As  before 


(n  +  p)  3  yb  and  Agb  3  0  hence  Agn  3  -AgP  approximately.  This  represents  four 
equations  in  five  unknowns.  However  it  is  known  that  nQ  3  0  so  one  of  the 
unknowns  is  eliminated  and  we  can  discard  the  first  column  of  AQ  to  produce 

-1  0  0  0 

O  42  -43  0  0 

B0  * 

U  0  6-70 

0  0  8  -1 


and 

_  -T 

8gn  -  “AqP  »  n  -  (n.j 

- 1 

Now  we  have  four  equations  in  four  unknowns  and  BQ  exists,  but  unfortunately 
- 1 

Bq  is  not  an  integer  matrix.  Again  start  factoring. 

1  0  0  0  1  0  0 

-42  1  43  0  0  301  0 

0  0  1  0  0  0  1 

0  0  0  1  0  0  0 


Next  left  multiply  both  sides  of  Bgn  =  -AqP  by 
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42/301 

1/301 
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[  Again  force  the  lefthand  side  to  have  integer  components. 

i  i 


Now  the  4x4  matrix  on  the  left  side  does  have  an  integer  inverse  and  the 


S  solution  to  this  system  of  equations  is 

f  :  ' 


.33 


